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1. Introduction and Summary 

Strongly interacting qnantum field theories when coupled to finite magnetic field exhibit 
a host of interesting phenomena^. In QCD such effects include the modification of the 
phase diagram of QCD at finite B [2], the chiral magnetic effect [3,4], and the magnetic 
catalysis [5], or de-catalysis [2] among others. These problems are much beyond mere 
academic interest as strong magnetic flelds B are produced in all non-central heavy ion 
collisions (i.e. those with nonzero impact parameter b) by the charged “spectators” (i.e. the 
nucleons from the incident nuclei that “miss”, flying past each other rather than colliding). 
Indeed, estimates obtained via application of the Biot-Savart law to heavy ion collisions 

^We refer to [1] for an extensive study of the various phenomena mostly for QCD. 
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with 6 = 4 fm yield e\B\/m^ ~ 1-3 about 0.1-0.2 fm/c after a RHIC collision with ^/s = 
200 AGeV and e\B\/m^ ~ 10-15 at some even earlier time after an LHC collision with 
^ = 2.76 ATeV [3,6-12]. 

We study the dependence of thermodynamics and the CP-odd transport in QCD in the 
deconfined phase at finite temperature and magnetic field in the limit of large QCD coupling 
constant and large number of colors W- In this limit, the gauge-gravity duality [13-15] 
allows one to study the theory by mapping it to a gravitational theory in one-higher dimen¬ 
sion. In particular we employ the bottom-up holographic model put forward in [16-19] to 
model the glue sector of large-A"c QCD at large’t Hooft coupling. However, this is not suf¬ 
ficient to study the effects of the magnetic field on the system as B couples the quark-gluon 
plasma through the quarks that constitute the flavor sector in the fundamental represen¬ 
tation of the gauge group SU{Nc), not the glue sector that is in the adjoint representation. 
The dynamics of flavor degrees of freedom are introduced by considering Nj pairs of flavor 
branes and anti-branes in the glue background, [20,21]. Then, non-trivial dependence on 
B of any quantity in the large Nc limit will only be visible if we also consider large number 
of flavors Nf and keep the ratio fixed: 


Mf —)> oo. 


Nf 

—A = X = const. 
Nc 


( 1 . 1 ) 


In the gravity dual, this limit requires back reacting the flavor branes on the background 
solution [22-27] and the solution now acquires a non-trivial dependence on the ratio x 
above. Dependence on any quantity on B will arise in the form xB. In this paper we 
consider finite but small x. This approximation, even tough makes our study slightly 
unrealistic (as x = 1 for QCD with 3 flavors), simplifies the calculations drastically as 
explained below. Therefore, in this work, we shall mostly confine ourselves in the qualitative 
features of the B-dependence in the system. Previous holographic studies of the A/" = 4 
Super Yang-Mills thermodynamics in the presence of magnetic field include [28]. The 
sphaleron rate in AA = 4 SYM, for finite B was calculated in [29]. 

We study two separate effects of the magnetic field in this paper. Firstly we ask 
how the entropy density S of the thermal state depends on B. We find that the ratio 
S{B Q)/S{B = 0) increases with B at any fixed temperature above Tc. Here Tc is 
the deconfinement temperature. We also find that the rate of increase becomes more 
substantial at lower temperatures T > Tc and the dependence oi S{B 0)/S{B = 0) on 
B becomes milder as T is raised. This finding is summarized in figures 1 and 2. All of 
this is in non-trivial qualitative agreement with the recent lattice studies [2]. Secondly, we 
consider anomalous transport properties of QCD at hnite B and T. In particular we study 
how the sphaleron decay rate (sometimes called the Chern-Simons decay rate) depends on 
these variables. 

In QCD at finite temperature there exist sizable effects, e.g. the sphaleron decays [30] 
that generate a non-trivial expectation value (Tr Gf\G), which in turn generates an effective 
chiral chemical potential /xs [3,4] for the quarks. This effective thermodynamic variable is 
an important ingredient in the recently discovered Chiral Magnetic Effect [3,4]. In short, 
the CME refers to generation of a macroscopic electric current in the presence of an external 
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Figure 1: [Color online] Top: the entropy density S, as a function of eB with IHQCD 
potentials. Down: change in entropy density divided by vs. T/Tc for different values of 


eB. 


magnetic field B in gauge theories with chiral fermions as a result of the chiral anomaly. 
A number of independent derivations [4,31] reveal that the generated electric current is of 
her form 

J = o-cmeB, ( 1 . 2 ) 

where the so-called chiral magnetic conductivity ctcme is of the form 

e2 

(7CME = M5 • (1-3) 

Therefore, it is essential to determine the dependence of /rs on B and T in order to assess 
the importance of this phenomenon. On the other hand, the most effective mechanism 
that generates in QCD are sphaleron decays [30]. Therefore the question translates 
into a calculation of the sphaleron decay rate in QCD at finite B and T. This rate in any 
Quantum Field Theory is captured by the Wightman two point function of the topological 
charge, 

^ ^ ^ (1.4) 
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Figure 2: [Color online] Tqs eB with IHQCD potential for different temperatures. 


With Cl = 1 and C 4 = 1. 

where In a state invariant under translations in space and time, the rate 

of change of A^cs per unit volume V per unit time t is called the Sphaleron decay rate, 
denoted Fes, 



(1.5) 


where the subscript W denotes the Wightman function. Therefore, in order to study the 
magnitude of the chiral magnetic current, one should study this Wightman correlator of 
the topological charge. Holographic calculations of rc 5 in different holographic models are 
presented in [32,33,36], In the regime of strong interactions, this quantity can be obtained 
by means of the gauge-gravity correspondence by studying the propagation of a bulk axion 
field in the 5 dimensional gravitational background. Our findings are shown in figure 2. We 
find that, similar to the entropy above, the sphaleron decay rate increases as a function of 
B and the rate of increase is more substantial at low temperatures. Thus, we conclude that 
presence of the magnetic held in QGP magnihes the chiral magnetic effect, in an indirect 
fashion, by increasing /rs in (1.3), in addition to the direct dependence on B in (1.2). 

The rest of the paper is organized as follows. In the next section we introduce our 
holographic model that corresponds to QCD with havors in the Veneziano limit (1.1) and 
present the equations of the motion for backgrounds at hnite T and B. In this section we 
also introduce the CP-odd part of the gravity action, that contains the axion held, that 
will be used in calculating the Sphaleron decay rate. In section 3 we present the calculation 
of the background that we obtain by solving the background equations numerically for a 
specihe choice of the model. In the same section we present our hndings for the B and 
T dependence of thermodynamic variables, in particular they entropy density S{B,T). 
Section 4 is devoted to the sphaleron decay rate, that we calculate by solving the axion 
huctuation equation numerically, on top of the numerical background found in section 3. 
We present our results for dependence of the sphaleron decay rate on B and T in this 
section. Section 5 is a discussion and an outlook on our research. Two appendices present 
details of our calculations. 

Note added: Another paper [34] that studies the dependence of thermodynamics on 
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magnetic fields appeared on the arXiv yesterday while this work ready to post. While we 
have some overlap with this paper in our results concerning entropy as a function of B, 
in general we consider different holographic models and the focus of the two papers are 
somewhat different. 

2. Holographic QCD with Flavors 

Improved Holographic QCD, [16-19] , is a string theory inspired bottom-up model of large 
Nc 4-dimensional pure Yang-Mills at strong coupling, which is in remarkable agreement 
with low energy QCD phenomenology. The model was generalized in order to include 
mesonic physics in the Veneziano limit, where Nf —)• oo, Nc ^ oo and Nf/Nc = finite, 
A = g^Nc = hxed, [22,23,25-27] . The full action for the Veneziano QCD model (V-QCD) 
can be written as 


S = Sg + Sf + Sa ( 2 . 1 ) 

where S'g, 5/, and Sa are the actions for the glue, the flavor and the CP-odd sectors 
respectively. As discussed in [16] and [ 22 ], only the first two terms contribute to the 
vacuum structure of the theory if the phases of the quark mass matrix and the 9 angle 
vanish, hence the ground state is CP-even. In this work, we study the theory above the 
conhnement-deconfinement phase transition in the presence of an external magnetic field 
which descends from the flavor degrees of freedom. However, Sa is important for the 
analysis of the CP-odd excitations of the model. 

2.1 The glue sector 

The holographic glue action is a two-derivative gravity-dilaton action with fields that cor¬ 
respond to the lowest dimension operators of the boundary field theory. The metric is dual 
to the energy momentum tensor of the theory and the dilaton corresponds to the TrT^ 
operator. The action was introduced in [16] and reads 

Sg = I + P3(A)) . (2.2) 

Here A = is the exponential of the dilaton field and its boundary value is identified with 
the holographic’t Hooft coupling. Both the metric and the dilaton field are non-trivial in 
the ground state^. 

2.2 The flavor sector 

The flavor action is the generalized Sen’s action which was first used in holographic QCD 
in [ 20 ], in the probe limit, and was incorporated in the study of backreacting flavors in [ 22 ]. 
The form of the action is 

^Our choice for the dilation potential as well as the other potentials in the action is presented in section 
3.1. 
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Sf = -^M^NcTr j d^xdr (Vf{X, T^T)^- det Al + Vf{X, TT^)^- det Ar) , (2.3) 

The quantities inside the square roots are defined as 

{DmT)\DnT) + {DnT)\DmT)\ , 


^LMN = 9MN + w{X, T)F^lf H-^^ 


A-rmn = Qmn + w{X, TjF^j^ H--— 


with the covariant derivative 


{DmT){DnT)^ + {DNT)iDMT)A , (2.4) 


DmT — OmT + iT— iA^T. 


(2.5) 


The gauge fields A^, Ar, and the complex scalar T transform under the U{Nf)^ x U{Nf)ii 
flavor symmetry and they are dual to the left/right axial current and the quark mass 
operator, respectively. One can define the vector and axial vector combination of the 
gauge fields 


Vm = 


aL aR 


aL _ aR 

A _ 

, Am — • 


2 ’ 2 
We also define the ratio of the number of flavors over the number of colors 

Nf 

The form of the tachyon potential is generally 

Vf{X,TT^) = V/o(A)e-“(^)^^^ . 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


The Plank mass, M, which appears as an overall factor in front of both Sg and 5/, is 
fixed by requiring the pressure of the system to approach the large temperature limit of a 
system free non-interacting fermions and bosons. This fixes (Mi)^ = (l-|-7x/4)/457r^, [23], 
where i is the AdS radius. The coupling functions k{X,T) and w{X,T) are allowed in 
general to depend on T, through such combinations that the expressions (2.4) transform 
covariantly under flavor symmetry. However, following string theory intuition we take 
them independent of r. The potentials Vfo{X), a(A), k{X), and w{X) are constrained by IR 
properties of QCD like chiral symmetry breaking and meson spectra as it is studied in [25]. 


2.3 The CP-odd sector 

The CP-odd part of the action has been studied in detail in the probe limit in [20]. Its 
main features are that it arises from the Wess-Zumino action that couples the Ramond- 
Ramond forms with the gauge fields on the branes and it is such that it reproduces the 
U{1)a anomaly of the held theory. For hnite x, the action was introduced in [25], and it 
couples the axion held, dual to TrF A F, from the closed string sector to the tachyon’s 
phase (T = re^^lNf) and the U{1)a axial-vector held, in such a way that reproduces the 
correct U{\)a axial anomaly of the boundary held theory 
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Sa^_M^ I [da-x{2Va{\T)A-^dVa{\T)f , 

where ^ is the phase of the tachyon field. The action is normalized so that a is dual to 
9/Nc with 6 being the standard 6*-angle of QCD. The potential Z{X) has been studied in 
the X = 0 case extensively and is constrained by the topological susceptibility of QCD and 
the 0"' glueball spectrum [17]. In this work, we are not interested in the exact form of the 
potential Va- We notice though that it should satisfy Va{X,T = 0) = 1, in order to have 
the correct U(1 )a anomaly. A detailed analysis of this action in the zero temperature case 
is under preparation. 

2.4 Background at finite magnetic field and temperature 

The non-trivial bulk fields on the ground state include the metric, the dilaton, the tachyon 
and a constant magnetic field. The non-zero background magnetic field arise from the 
flavor part of the action. Hence, it is naturally defined as the U{1) part of the vector field 

Cm= (0,-yB,^H,0,0) , Am = 0. (2.10) 

The non-trivial magnetic field breaks the SO{3) rotation symmetry of the vacuum state to 
50(2) symmetry of the directions transverse to the magnetic field plane, — x^. Hence, 
the Ansatz for the background metric and the dilaton should be, 

ds^ = - f{r)dt^ + dxl + dxl + , A = A(r), (2.11) 

The UV boundary lies at r = 0 (and A —)■ oo). In the UV, the AdS coordinate, r, 
is identified roughly as the inverse of the energy scale of the dual field theory. At zero 
temperature, /(r) = 1 and the solution has AdS^ asymptotics near the boundary with 
logarithmic corrections. In the IR region of the space-time the solution asymptotes to a 
qualitatively similar solution for A{r) and A(r) as the B = 0 case, while the function W{r) 
will be proportional to B. At finite temperature, at least two solutions should exist. One 
is the similar to the zero-temperature solution but with periodic time coordinate and a 
black hole solution. At some certain temperature, confinement-deconfinement transition 
is expected to happen, similarly to AA = 4 case, [28]. Above the deconfinement transition 
the dominant solution is a black hole metric, so f{r) is non-trivial and satisfies f{rh) = 0, 
where r = r/i is the position of the black hole horizon. The black hole temperature is 
T = /3“^ = j/'(rh)j/47r. The boundary asymptotics should match the T = 0 solution. 

Moreover, in the general case of a boundary field theory with Nf light quarks, we 
consider a background tachyon field of the form 

r = r(r)Ijv,, (2.12) 

which corresponds to quarks with the same mass. Since, r(r) is dual to the quark mass 
operator, its profile signals the chiral symmetry breaking in the boundary field theory. 
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Background solutions with non trivial tachyon which is diverging in the the deep IR of 
space-time generically correspond to chirally broken state of the field theory. Bulk solutions 
with identically vanishing r(r) signal a chirally symmetric phase of the boundary theory. 
V-QCD, at R = 0, has been studied both at zero and non-zero temperatures and it has been 
found that the theory exhibits a chiral transition at finite temperature, depending on x, 
above which the chiral symmetry is restored, hence r(r) = 0 and the theory is deconhned. 
In the present work, we analyze the theory in the chirally symmetric phase and non-trivial 
B. A full study of the model would require to find all the different bulk solutions, with the 
same near boundary conditions, for each B. By comparing their free energies one can decide 
which is the dominant solution. In the present work, we assume that the magnetic field 
does not change the phase structure dramatically, so at high temperatures the dominant 
phase is deconhned and chirally symmetric^. The CP-odd helds are set to zero at the 
ground state, a = ^ = Am = 0. The vacuum action then reads 


S = / d^x 




V 3 A2 

: V/(A, r) ^J-det{g^y + w{\) + k{X, t) d^r d^r) 

The Einstein equations of motion read 


— X 


(2.13) 


Rfiu 


— X 


2' 

VfiX,T) 


4 dfj,XduX 


3 A2 


2{dXf 1 ,, ,,, 

dD . 

= 0 , 


^/D dgt^^ 


(2.14) 


where D = det{5^ -|- w{X) V^x -|- k(A, T)g^’^duT d\T). The full set of equations of motion 
for non-zero r(r) is presented in Appendix A. We now consider the deconhned and chirally 
symmetric phase where the fermion condensate is zero and we have a background magnetic 
held. This magnetic held comes from the vector held on the havor branes and is described 
by Sen’s action (2.3). In case of small magnetic held the DBI can be expanded to quadratic 
order to the Maxwell action 


Sf = —M^ NcTr f d^x drVf{X, det((5(l -|- w{X, TYg^PVpy) 

r / lA 12 X (2-15) 

= -m3 W Nf J d^x dr Vf{X, r) ^ ( 1 + j . 

The gauge held part of the action reads 

Sb = -m3 W Nf J d^x dr V > (2-16) 

®This can be further justified since the condensate in QCD decreases for increasing magnetic field at 
temperature close to Tc, a phenomenon called inverse magnetic catalysis, [35]. 
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where Vb(A) = Vf{X)w{X)'^. The potentials are taken independent on the tachyon since it 
is neglected in the current solution. We notice that for small x the contribution of the 
magnetic field in the background solution is small. The Einstein equations now read 


Rfiu 1 2 ^2 


4 d^XdyX 2 {dX) 


^2 9^lv + 


— X- 


VbiX) 


y^^Vup - ] = 0. 


(2.17) 


where now the dilaton potential is replaced by Veff = Vg — xVjq . The equations of motion 
for the metric ansatz functions are 


3A"{r) + l^ + 3A'^ + {3A' - W') 


2 „-2A 


o2A 


3A^ ' ' '.'2/ ' 4/(r) 


W'f 


/w/ xVb{X)B e 


2 ^-2A 


W" + + W'2 + 3MW' - ^ ^ , 

/ 2/(r) 

/" + (SAl' + W') f -X VbiX) e-2^ = 0 . 

We can integrate the equation of /(r) and find 


= 0 , 


fir) = e 


2W 


1-Ci 


e-^^-3Wdr' 


(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 


where we considered the boundary conditions, /(O) = 1 and kE(0) = 0. Ci is determined 
by requiring regularity on the black hole horizon 

1 


Cl = 




Using 2.21, this integration constant can be related to the enthalpy as 

Cl 


5 = 


The first order constraint equation reads 


IbvrGs T 


( 2 . 22 ) 


(2.23) 


X 14(A) ^2 e-2^ 
4/ 


A f'iP) a a>2 ^ . e 


2A 




= 0 


(2.24) 


The dilaton equation of motion is not independent but can be derived from combining the 
above equations 


} —\/ \2 4” ( {r) + lU'(r) + ' I . ) ! + - ^ \ - d\VeffiX) 

A(r) A(r)2 V f{r) J A(r) 8 /(r) 


+ (^324'(r) + W'{r) + 


f'{r)\ A'(r) 3 A(r) 


3B^e 2"^U)A(r) 


9AUb(A) = 0. 


16/(r) 

The equations of motion enjoy the following scaling symmetries 


(2.25) 
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• rescaling of /(r), / AA - ^log Cf , B ^ 

• rescaling of the AdS coordinate, r ^ Ar , A ^ A — log A , i? —^ , 

• rescaling of the W{r), W W + cw- 

2.5 UV asymptotics 

The UV asymptotics of the above equations that we are interested in are AdS with loga¬ 
rithmic corrections [16,17]. The magnetic field influences the UV expansions of the fields 
at O(r^). In case of a very strong magnetic field and a potentials 14// , 14 of the form 


14 // ~ + +^2^^ + ...) , 14 ('^) — + UftiA + ...) (2.26) 



where Ct is a temperature dependent coefficient. 

2.6 IR asymptotics 

Requiring regularity on the horizon, the near horizon asymptotics of the system are of the 
form 


A = Ah + A!j^e + , W = Wh + W'^^e + ... 

A = A/i + e -|-... , / = /^e -|- + ... , 


(2.28) 


where e is the distance from the horizon, e = — r. Expanding the equations of motion, 

(2.20, 2.24, 2.25) we find the coefficients 


r 2 ^X. i '^Vb [Xh) 

A. 



fh = -e^^^VeffiXh) + h-^^'^B^xVbiXh ). 


(2.29) 


3. Numerical solution 


3.1 Choice of the model 

In Improved Holographic QCD, the asymptotics of the potential 14(A)^are fixed for small 
A (UV of the field theory) to match the perturbative large-Vc /3-function, which gives the 













dual field theory asymptotic freedom. In large A (IR of the field theory), the potential must 
have a form such that the dual field theory is confining and has a linear gapped glueball 
spectrum [16-18]. In addition, as r —)■ 0, we want metric to approach AdS, ^ and 

A —)• —l/log(r) to mimic the perturbative running of large-A'c YM coupling. A form with 
the correct asymptotics that we use in our numerical calculations is 

R(A) = ^(^l + VoX + RiA^/yiog (1 + R2A4/3 + . (3.1) 

We can constrain Z{X), the dilaton-dependent normalization of the axion’s kinetic term, 
as done in [36]. We will use the following forms for Z 2 with correct asymptotics, 

Z(A) = Zo (1 + C4A^) (3.2) 

Z{X) = Zq (1 -|- ciA + C4A^) , (3-3) 

with coefficients ci and C 4 restricted by lattice calculations of axial glueball masses [36] as 

0 < Cl < 5, 0.06 < C4 < 50. (3.4) 


The potentials appearing in the flavor action, (2.3), can be constrained by looking at 
their asymptotics and comparing to lattice and perturbative results as done in [22,23,25]. 
Their UV asymptotics are chosen to match the perturbative anomalous dimension of the 
quark mass operator. Their IR asymptotics were fixed by QCD features of the flavor sector 
as chiral symmetry and the meson spectra, [25]. Here we simply present the form of R/o(A) 
and w{X), that we use in our calculation 


w{X) 

VfoiX) 


(1 + 3|iA)^/3’ 

Wo (1 + WiA + WsA^) , 


(3.5) 

(3.6) 


with coefficients 


ai 


ITi 


115 - 16x 
2167r2 ’ 

24-h (11 - 2x)Wo 
277r2Wo 


Wo = 3/11 


W 2 


24(857 - 46x) -b (4619 - 1714x + 92x'^)Wo 
466567r'^lTo 


(3.7) 

(3.8) 


3.2 Numerical technique 

In order to define a physical, finite action for a non-compact geometry, we must choose a 
reference background with the same asymptotics for the metric as well as the dilaton. 

We solve Eq. (2.17) using the following numerical method. First, we note that the free 
integration constants of this solution are //, A^, Wh-, 4>h = logA/j (letting the subscript h 
denote the function at r/^), as all other integration constants are fixed in terms of these 
by the Einstein’s equations and the condition of horizon regularity, (2.29). The only other 
free parameters are the location of the horizon, and B, the magnetic field strength. We 
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choose the reference background to be the B = 0 solution. Then, at some UV cutoff Tc, for 
each value of 5 7 ^ 0 , we match the geometry of the solutions at the cutoff by demanding 


V/(u)/3 = V/o(u)/3o, 

II 

0 

(3.9) 


V 0 I 2 = V 0 I 2 

(3.10) 

A(rc)+H^(r-c)voli = 

1 

II 

(3.11) 

(t>{rc) = 4>o{rc)- 


(3.12) 


where the subscript 0 means B = 0. The left column contains the requirements for matching 
the intrinsic Euclidean geometry of the two solutions and the right column contains the 
extra conditions we add for convenience. For example, it is convenient to see how a solution 
changes with B while keeping temperature T = 1/(3 = fixed. 

Our numerical method is then a shooting method, that is, we vary the free horizon 
quantities Ah, Wh, (l>h, fh until the constraints at the cutoff are satisfied. Now, we can see 
how the horizon quantities Ah, Wh, <i>h, fh change as a function of B. 


3.3 Solutions 


We can use the UV matching procedure described above with a choice of x = 1/10 to solve 
for Ah, Wh, 4>h, and Vh- In the figures below, we present the dependence on the physical 
magnetic field eBphys where e is the elementary charge and Bphys, the physical magnetic 
field of our dual field theory, is obtained by 

Bphys — (3.13) 

where £ can be found in physical units by taking (. k, °' 247 MeV^ ’ 

Once the background solution is found, it is a trivial matter to calculate the entropy 
of the background using the Bekenstein-Hawking entropy formula 


Area/i Q^M'^h)+W(rh) 

4G5 ^ 


(3.14) 


where Area/^ denotes the area of the horizon and G 5 is Newton’s constant, related to the 
5-dimensional Planck mass by = l/( 167 rG 5 V^). This can be converted into a function 
of T and B as the location of the horizon Vh can be numerically obtained for varying values 
of B and T. An important point here is that the dependence on B and T is twofold: in 
addition to the explicit dependence on Vh, hence T and B, in (4), the metric functions 
A and W also change as we change B and T resulting in a hidden dependence on these 
variables. 

The plots^ in Fig. 3 show the change in entropy as a function of eBphys for different 
temperatures. In addition we plot in Fig. 4 the entropy’s dependence on temperature, for 
different values of eB. 

It is tempting to compare these results with the recent lattice studies, see for example 
figure 10 in [37], that we reproduce here in figure 5. Comparison of figures 5 and 4 show 

®Tc in the figures 3 and 4 below refer to the deconfinement temperature at B = 0. 
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Figure 3: [Color online] The entropy density S, as a function of eB with IHQCD potentials 
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Figure 4: Change in entropy density divided by vs. T/Tc for different values of eB. 


very good qualitative agreement. To see that the agreement is better than it appears in 
these figures one has to consider only the part T > Tc of figure 5. The tails on the left 
of the peaks at around Tc ~ 150 MeV in this figure are absent in the large Nc limit that 
we consider here. One can even hope for quantitative agreement and the reason that our 
results do not quantitatively agree with that of the lattice studies should be because we 
consider a model where the ratio of flavors to color x = 0.1, whereas in [37] this ratio is 
X = 1. We expect the agreement become better when we consider larger values of x. 

3.4 Free energy and phase transition temperature 

Another important thermodynamic quantity is the free energy. In the case where the en¬ 
tropy only depends on the temperature, the free energy can be obtained from the entropy 
that we found above simply by integrating the first law in T, [18]. On the other hand, this 
cannot be done here as 5 is a function of both B and T. Then the only way to calculate 
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150 200 250 

T (MeV) 


Figure 5: [Color online] Change in entropy density divided by in terms of T/Tc for 
different values of eB. It is noticed that it has a maximum at Tc- The plot was produced 
using Lattice field theory in [37]. 


the free energy from the holographic dual theory is from the gravity action (including the 
Gibbons-Hawking term and the counter-terms) evaluated on the background solution [18]. 
A similar method is to consider the difference of the on-shell actions evaluated on the black- 
hole and the thermal gas solutions, that is a solution with the same boundary asymptotics 
but no horizon. As explained in [38], (see also [18,19]) the black-hole solution corresponds 
to the deconfined plasma phase and the thermal gas corresponds to the low tempera¬ 
ture, confined phase. The curve in the (i?, T) space where this difference vanishes then 
determines the phase boundary in the system on which the confinement/deconfinement 
transition takes place. This curve Tc = Tc{B) can be calculated in the aforementioned 
manner. 

We can easily obtain the thermal gas solution that obeys the same boundary asymp¬ 
totics as the black-hole solution above by employing the following trick. The thermal gas 
has no horizon hence vanishing entropy. Using equation (2.23) then this means setting 
the integration constant Ci in the solution (2.21) to zero. Hence the thermal gas solution 
obeys 

/(r) = exp(2kF(r)). (3.15) 

The other metric functions can then be determined by numerical integration. Once both the 
black-hole and the thermal gas solutions are obtained, one should calculate the difference of 
the on-shell actions to determine the phase diagram. We shall not carry out this calculation 
in this paper for two reasons: firstly there are various numerical difficulties which render 
the calculation very tricky and one needs new techniques to maintain numerical efficiency 
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of the difference of the actions. Secondly, we consider small values of x in this paper, hence 
the difference Tc{B) —Tc{0) = 0(x) that is small. This will make a difference 0{x‘^) in the 
figures 3 and 4 above, which we can safely neglect. 


4. CP-odd fluctuations at flnite x and magnetic fleld B. 


We now analyze the coupling of CP-odd excitations of the model and calculate the Chern- 
Simons diffusion rate at finite B. As explained in section 2.3, the pseudo scalar axion a is 
dual to the topological charge operator given by 


Recall that the Chern-Simons diffusion rate is given by 


rc 5 = lim lim 

V—^-CXD t^CO 


Vt 


(4.1) 


(4.2) 


This can be rewritten as 

'^CS = J {q{x>^)q{0))w, (4.3) 

where {)w denotes the Wightman correlator, [30]. 

In Fourier space, let Gw{^-,k) and GR{uj,k) denote the Wightman and retarded 
Green’s functions of the topological charge operator q, defined in Minkowski spacetime 
as 


Gnik) = -ij d^xe-^^--e{t){[q{x),qm) (4.4) 

Gw{k) = \ j d‘^xe~^^'^{q{x)q{0) - q{0)q{x)) (4.5) 

where k = (a;, k). Then, the fluctuation-dissipation theorem relates these as 

Gw{uJ,k) =-coth(^-^^lmGR{u},k). (4.6) 

for CO <^T. Taking the zero momentum and small frequency limit of this, we can write the 
Chern-Simons diffusion rate as 

2T ~ ^ 

Tcs = lim — Im Gr{uj, k = 0), (4.7) 

ui 


where we now consider Gr{u!, k) to be the retarded Green’s function of the topological 
charge operator q{x^), dual to the axion. To compute the above transport coefficient in 
our model at finite B, we study the quadratic excitations of the axion field which are 
coupled to the U{1)a axial-vector current and the phase of the tachyon. The interplay 
of the gluon topological correlator and the axial current has been studied in the probe 
limit in [39-41]. In the DBI, (2.3), the coupling of the axial-vector to the phase of the 
tachyon is proportional to the tachyon background solution which is taken zero, hence 
this is neglected. At r(r) = 0, the potential appearing in the CP-odd action, (2.9), is 
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V"a(A,r = 0) = 1, so the axion coupling to the phase of the tachyon is also zero since it is 
proportional to the derivative of 14- Expanding the actions (2.3, 2.9) up to second order 
in the fluctuations we have, 


5i = I d^xdrVbiX) V-det(g) AstGs™ A^m 


1 


+ Ga’^^ Ast Ga^^Anm - - {Ga^^Amn)'^ 

^2 = — J d'^xdry/^ Z{X) (yduoi — 2x AM){dis[Oi — 2xVa . 

Si conies from the expansion of the DBI and the metric G is defined as 

/-e2^W/(r) 0 0 0 0 \ 

0 w{X,t)B 0 0 

0 -w{X,t)B 0 0 

0 0 0 f,2A{r)+2W{r) q 


Gmn = 


0 


0 


0 


0 


/4) ’ / 


(4.8) 


(4.9) 


It is noticed that Gmn is split in its symmetric and antisymmetric parts Gmn = Gsmn ~X 
Gamn- The equations of motion read to leading order in the magnetic field 


Qm 


On 


Z{X)y/^g^^{dMa{r,x)-2xAM{r,x)) =0 
Vb{X) i/-det{G) ( Gs^^Gs'^^Ast + Ga^^Ga^^ Ast 

+ 2Z{X)^f^g^^ dNa{r,x) -2xVaAN{r,x) 


(4.10) 


+ ^Ga^^Ga^^Ast 


= 0 


In order to calculate the Chern-Simons diffusion rate we now look at the equation of 
motion for the axion fluctuation a in the small lo limit. Following [32], we solve this 
equation with a Dirichlet boundary condition at the asymptotically AdS boundary (at 
r = 0) and an in-going wave boundary condition at the horizon r/j. Considering the case 
where /c^ = (—w, 0, 0, 0), the solution for the axion takes the form 


a{r,t) = j *^*a(r,w)a(w). 


where a{oj) is fixed by the Dirichlet boundary condition, 

linia(r,t)= [ ^e"*‘^*a(a;), 
r —^0 / ZTT 


(4.11) 


(4.12) 


and a(r,uj) is the solution to Eq. (4.10) in momentum space. Similar expansions hold for 
Am- Taking the spatial components of the axial-vector field to be zero, the fluctuation 
equations become 
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dr 


dr 


Z{\) e 


3A+W 


fdra 


ioj Vfc(A) dr At 


ioj V'fe(A) dr At - 2Z(\) f drO = 0 


^3A+W 

+ Z (A) — j — — 2iiox At) = 0 

„3A+W 

+ 2Z{X) — y — (w^a — 2iujx At) = 0 

3A+W 


(4.13) 


In the present work, we restrict ourselves in the calculation of the Chern-Simons diffusion 
rate, which is the transport coefficient of q{x), hence we are interested in solving the 


dr 


Z{\)e^^+^ fdro 


p3A+W 

+ Z(A) — y — Lo^a = 0 . 


The on-shell action of the axion reduces to the boundary term 


iJ. 


on-shell 


= f — a{—u}) T{r,uj) a{u!) 


2tx 


rh 


where 


m3 


T{r, to) = —-^) Z(Kr)) dra(r, u). 

The retarded Green’s function as prescribed by [32], is given by 


(4.14) 


(4.15) 


(4.16) 


Gr(oj) = —2 lim J^(r, co). 

r—tO 


(4.17) 


In order to calculate rc' 5 , we need to solve Eq. (4.14) for a with small oj. We do this 
using near-horizon matching following [36]. In this method, we first solve Eq. (4.14) with 
a; = 0 and expand the solution near the horizon. Then we reverse the order, solving in the 
near-horizon region and then expanding in small ca. Einally, we match these to solutions 
to obtain T. 

When cu = 0, the solution of Eq. (4.14) is 


5a 


pr 

Cl+ C 2 

Jo 


dr' 

Z(A(r'))e3^h')+'^h')/(r')’ 


(4.18) 


with Cl and C 2 constants. The integral term in Eq. (4.18) diverges as r —)■ r/^, since 
fii^h) = 0. Therefore, a normalizable solution must have C 2 = 0 when w = 0 and C 2 oc a; 
when u! is small. This solution Eq. (4.18) gives for Eq. (4.16), 


lim /"(r, uj) 
r—tO 



uj<^T, 


(4.19) 


In order to have unit normalization at the asymptotically AdS boundary, we choose Ci = 1. 
Now, expanding Eq. (4.18) near the horizon, we find the solution 


5a = Ci + 


C 2 


Z{\h)e^A^^A+W{rt.) f^ru) 


log(rh - r) 0{rh - r), 


(4.20) 
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Figure 6 : [Color online] Fc^ vs. eB with IHQCD potential for different temperatures. 
With Cl = 1 and C 4 = 1. 


where we used \f'{r}^)\ = dvrT. Next we do the reversed order of operations. First, we 
expand Eq. (4.14) near the horizon and we find the solution 

iu) _ iu) 

a = C+(r/i — r)47rT + — r) 4,nT ^ (4.21) 


with coefficients C+, C- depending on oj but not r. In order to have an in-going wave 
condition at the horizon, we set (7+ = 0. Expanding the solution for small cj, we get 

^ = ^-~ -r) + 0 . (4.22) 

By matching the constant and logarithm terms in Eqs. (4.20) and (4.22), we find 

Cl = C_ = 1, C 2 = (4.23) 


Using Eqs. (4.17) and (4.19), we obtain 

Gr{uj, k = 0) = -iwM 3 Z 2 (A/,)e 3 ^('''‘)+^('’'^), a; < T. (4.24) 


And so our final result for Fes is 


^cs 


1 T 
i^SvrCs 




1 sT 


K^Z{Xh), 


(4.25) 


using that the entropy density is s = and = l/( 167 rG 5 A"^). 

The dependence of rc >5 on temperature in Fig. 6 is in good agreement with the previous 
result in [36], showing a peak at the devonfinement transition. We notice that rc 5 is 
quadratic in B for weak magnetic field and becomes linear for strong magnetic field. 
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5. Discussion and Outlook 


In this paper we studied the effect of finite magnetic field on the quark-gluon plasma in the 
Veneziano limit (1.1) at strong coupling, using a bottom-up holographic model for QCD. 
The Veneziano limit is crucial because the magnetic field couples to the plasma only through 
the flavor sector and in the large color limit, this coupling effectively vanishes unless one 
takes the number of flavors to infinity as well, keeping the ratio fixed. We found that the 
thermodynamic variables exhibit non-trivial dependence on B. In particular the entropy 
density increases with B and the rate of increase is more pronounced at smaller values 
of T, close to the deconfinement temperature T^. This result is in very good qualitative 
agreement with the recent lattice studies [37]. We also studied the sphaleron decay rate 
and found that this rate also increases with B. This means that sphaleron decay processes 
are favored in the presence of magnetic field which in turn imply a larger value for the 
effective axial chemical potential in (1.3). Therefore the chiral magnetic effect is more 
probable in a strongly interacting plasma with magnetic field. 

There are several future directions to explore. First of all, we simplified our calculations 
by taking a small value of x, the ratio of flavors to colors. This allowed us to expand the DBI 
action on the flavor branes. It is straightforward but cumbersome to lift this approximation 
and study the system in a more realistic case of x ~ 1. We leave this study for the future. 

Secondly, we have not studied the dependence of the free energy on B in this paper. 
As explained at the end of section 3.4, this, and the phase diagram of the theory can 
be explored by studying the difference of the on-shell actions evaluated on the black-hole 
(deconfined phase) and the thermal gas (confined phase) solutions. We plan to study this 
in the near future. In fact, one can also add the baryon chemical potential fi to the phase 
space and explore the full phase diagram in the space (T, B, fi). 

Finally, we are currently studying the response of the axial charge to the fluctuations 
of the topological operator, which requires the study of the full, coupled system of the 
CP-odd excitations. Similar issues are addressed in the probe limit in [39-41], 
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APPENDIX 

A. Equations of motion for r 7 ^ 0 

Einstein equations are written in terms of the ansatz functions of Eqs. (2.10, 2.11, 2.12) 
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3A"M + + 3^- + (3A' - 1 ) - = 0. 


2/ 2Q(r)/(r) 


2/(r) 


W"-^ + W"- + 3A'W^' + (1 - 0^) = ». 

/ 2Q(r)/(r) 


(A.l) 


/'' + (3A' + H^')/'-^2iA,^(l-<3^)=0. (A,2) 

where we have defined the functions 

G(r) = + e“2^WK(A, T)f{r){drT{r)Y , Q(r) = y^l + r(;(A, r)2i?2g-4yl(r)^ (A.3) 


The first order constraint equation reads 

2 \'2 fit \ 2 A 

__ _ (^^A' + W') -GA'^ - 3A'W' + —V,{\) 

xVf{X,T)Qe'^^ _ 

2Gf 


The dilaton eom is 


(A.4) 


A"(r) X'{rf 




+ 


A(r) A(r)2 
3x B'^ G{r) A(r) Vf{X, t)w{X, t) 

8/(r)Q(r) 

3 X A(r) Q{r) V/(A, r) r'(r)^ 


^ ^ 3a:e2^WG(r)A(r)Q(r) 

dxw{X,T) -- dxVf{X,T) 


8f{r) 


16 G{r) 


dxK{X,T) = 0. 


(A.5) 


This is not independent but can be derived from combing Eqs. (A.2, A.4) The tachyon 
equation of motion is 


r"(r) - 


g2A(r) 


/(r)K(A 


^9.1ogy^(A,r)+ e-2^M/(r)K(A,T) (w'{r) + ^^ 


^ A'(r)aAlog(/t(A,T)y/(A,r)2) - ^ log w{X,t)^ T'{rf 


+ (^A\r) 


+ W'{r) + + X!{r)dx log(V/(A, r) k(A, r)) - 


Q{rf 


A'(r)(l-Q(r)2) 


<5(0^ 


(9a log w{X,t) ) r'(r) 


+ log .(A,r) + d, log»(A,r) = 0 


f{r)K{X,T)Q{ry 


(A.6) 


To find the equation of motion of the gauge field, we define the matrix in the square root 
of the flavor action, Eq.(2.13), as 
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^LMN = ^RMN = QmN = 


0 
0 
0 
0 


0 0 
g2A(r) w{X,t)B 


0 

0 


0 

0 


0 
0 
0 

^2A{r)+2W{r) 

0 


0 
0 
0 
0 

e2A(r) -^2 

7 

(A.7) 


The inverse matrix is 


/- 


gMN _ 


\ 


.-2A(r) 

fir) 

0 

0 

0 

0 


0 

e-2A(r) 

Q(r)2 

Q(r)2 

0 

0 


0 

w{\,T)Be-^^^r) 

Q{r)2 

e-2A(r) 

S(r)2 

0 

0 


Then, the eom of the magnetic field is 


0 

0 

0 

^-2A(r)-2W{r) 


0 
0 
0 
0 

e-^A(r)jg^ 

G(r)2 , 


(A.8) 


dM (yf{K t) w{X, t) 7 ^ Ga^') = 0 > (A-9) 

where Ga is the antisymmetric part of G- The above equation is automatically satisfied for 
the gauge field of the form (2.10). 


B. Fluctuation equations r 7 0 and B^O 

The quadratic action of those fluctuations reads 


Si = X j (fx drVf{X, T)^/-det{G) f k(A, r) r(r)2 Gs^^ {dui + 2 Am) {OnC + 2Am) 


(B.l) 


w{X,t) 


2 r 


MS A n TN A I A n TN a (p MN a \2 

Gs AstGs AnM + Ga fiSTGA Majm--^\Ga mmn) 


S2 = j ct^xdr^f^Z{X) {dna - x{2Va Am - CdMVa) {dNa - x{2Va An - CdNVa))g^^ 


where the ’’metric” G was defined in (A.7). It is noticed that Gmn is split in its symmetric 
and antisymmetric parts Gmn = Gsmn ^Gamn- The fully couple system of the equations 
of motion in the presence of external magnetic field then reads 
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(B.2) 


Qm 

Qm 


Z{\) V^g^^idna - x{2Va Am - CdMVa) 


= 0 


Vf{X, t) k{X, t) y^-det{g) {On^ + 2An) 


(B.3) 


-Z{X)dMVaV^g 


MN 


OmO. + x{2Va Am — idMVa 


= 0 


Bn 


V/(A, r) Y^-det(^) w{X, r)^ ( Gs^^Gs^^A st + Ast 


STr mn a 
+ 2^-4 Ga ^st 


- 2 Vf{X, r) V-det(g) k{X, t) Gs^^ {Bn^ + 2An) 


+2Z (A) Va g 


MN 


BnOL — x{2Va An — ^BN^a 


= 0 , 


(B.4) 
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